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NOMENCLATURE 


A area 

a speed of sound 

B jet layer parameter 

Cy mass flow coefficient 

d diameter 

F(C) function describing resonant properties of the generator 
f(M) function describing acoustical properties of nozzle 
f frequency 

fo diaphragm-reed friction coefficient 

go acceleration of gravity, sea-level value 

H geopotential altitude 

ho ram-air jet width 

I flow Inertia parameter 

i = 

K non-dimensior^al cavity volume 

k[) diaphragm-reed spring constant 

A distance from ram-air jet to tube 

L tube length 

M( ) function describing resonant characteristics of the tube 
Mo molecular weight of air 

Ml generator flight Mach number 

mp diaphragm-reed effective mass 

m mass flow 

p pressure 

q ram-air jet velocity 

R reflection coefficient 



Nomenclature, (Continued) 


R* 

universal gas constant 

Rj 

radius of curvature, jet layer 
ind 

Re 

= — , Reynolds number 
Au 

ro 

effective earth radius 

T 

temperature 

t 

time 

U 

velocity 

V 

vol ume 

V 

velocity component across the flow direction 

W 

jet layer perimeter 

K 

downstream distance 


diaphragm-reed position 

1 

geometrical altitude 

Z 

= (o+iw) complex frequency 

a.P 

non-dimensional parameters 

Y 

specific heat ratio 

C 

non-dimensional complex frequency 

r\ 

diaphragm-reed non-dimensional parameter 


viscosity 
_ Jtij 

mj 

A 

ifiH 



tfiT 

a 

real part complex frequency 


Nomenclature, (Continued) 


diaphragm-reed non-dimensional parameters 
w imaginary part complex frequency 

SUBSCRIPTS AND SUPERSCRIPTS 

( )q corresponds to stagnation conditions 

( )x conditions in front of normal shock 

( )g conditions behind the normal shock 

(*) rate of flow 

( )' perturbation quantity 

( )„ free-stream flow conditions 

( )c corresponds to cavity conditions (See Figure 20) 

( )o relates to diaphragm- reed system 

( )e relates to conditions at the end of the tube 

( )h relates to flow conditions across exit holes (Figure 20) 

( )j corresponds to annular ram-air jet 

( )s relates to static flow conditions at the generator exit hoi 

( )x relates to resonance tube conditions 

( )w working equilibrium flow conditions 


ABSTRACT 


A combined experimental and analytical Invt^stlgatlon was performed to 
understand the various fluid processes associated with the conversion of flow 
energy Into electric power In a fluidic generator. Experiments were performed 
under flight-sImuiated laboratory conditions and results were compared with 
those obtained in the free-fllght conditions. From these measurements, It was 
concluded that the mean mass flow critically controlled the output of the 
fluidic generator. Cross-correlation of the outputs of transducer data 
indicated the presence of a standing wave In the tube; therefore the mechanism 
of oscillation was an acoustic resonance tube phenomenon and did not represent 
the Helmholtz type of oscillation phenomenon. 

On the basis of experimental results, a linearlied model was constructed 
coupling the flow behavior of the jet, the jet-layer, the tube, the cavity, 
and the holes of the fluidic generator. The analytical results also showed 
that the mode of the fluidic power generator Is an acoustical resonance 
phenomenon with the frequency of operation given by f« a/4L, where f Is the 
frequency of jet swallowing, a Is the average speed of sound in the tube, and 
L Is the length of the tube. Analytical results further indicated that 
oscillations In the fluidic generator are always damped and consequently there 
Is a forcing of the system in operation. 


X 


I, INTRODUCTION 


Periodic oscillations In jet flows Impinging on solid surfaces have been 
observed over a range of flow conditions and geometrical variations (Reference 
1). The oscillation phenomenon in such jet flows is very similar to that 
observed in axi symmetric cavity flows. Both experimental and an:.<>£lrfi.*l 
investigations of flows over cavities by V. Sarohia (References 2 to 6) have 
shown that the self-sustained oscillations are caused by the instability of 
the mean velocity profile and that they are amplified along the shear layer 
flow over the (javity. These previous investigations (Reference 5) further 
showed that the flow conditions downstream of the cavity influence the 
oscillation phenomenon as well as the shift cf flow oscillations from one mode 
to another. The particular mode of the oscillation Is determined by the phase 
of the propagating disturbance which attains the maximum value of the 
integrated amplification along the cavity shear layer. 

Experiments performed on axisymmetric cavity flow oscillations showed the 
oscillaticfi phenomenon depends on the external acoustic excitation (Reference 
7). The jet edge-tone configuration in the fluidic electric power generator, 
as sketched in Figure 1, is similar to one in the axisymmetric cavity flow. 
However, the jet edge-tone oscillations will be modified by the resonance 
tube, reed-diaphr«nm system, and cavity surrounding the fluidic generator* 

In flight, as indicated in Figure 1, air enters the fluidic power genera- 
tor through the port located in the nose and leaves through the exhaust ports 
spaced uni formally around the circumferences of the ogive. During the passage 
of the air through the generator periodic air flow oscillations are produced 
by the jet interaction with the resonance tube. These flow oscillations are 
transmitted through a diaphragm and rod to reed that switches magnetic flux 
within the coil, thereby inducing a voltage at the coil terminals. 
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Little is known, howeven, about the manner In which the amplitude and the 
frequency of these flow oscillations influence the conversion of flow energy 
into electric power in the system shown in Figures I and 7. The understanding 
of various flow interactions that occur in the fluidic power generator is very 
essential to the efficiency of conversion of flow energy into electric power. 
The investigation reported here was designed to advance the understanding of 
these unanswered fundamental questions. 

Laboratory measurements were made to relate to the free-f light 
performance of the fluidic generator. Based on the various underlying fluid 
processes associated with the production of oscillations in the generator, an 
analytic model which can predict the operation and the influence of 
configuration changes on generator performance was constructed. This model is 
the subject of the present report. The information obtained in this study has 
led to methods to improve the overall efficiency of the conversion of flow 
energy to electric power generation. 


II, EXPERIMENTAL ARRANGEMENTS AND MEASUREMENTS 


To enhance our understanding of the various flow processes responsible 
for the conversion of flow energy into electric power, two experimental setups 
were constructed. The bench test setup Is shown in Figure 2. This facility 
could supply compressed air at various stagnation temperatures up to 700®F. 

The air was heated by passing It through an electric heater with a temperature 
control Varlac as shown in Figure 2. This control led-temperature air was 
accelerated through a convergent nozzle which was flushed with the Inlet of 
the fluidic power generator. The total temperature and pressure In the 
settling chamber, along with static pressure at the Inlet were measured and 
utilized to compute the rate of mean mass flow entering the fluidic generator. 
The fluidic generator output, as Indicated In Figure 3, was monitored on 
the oscilloscope as well as on the rms meter. An environment chamber was also 
built to vary the ambient static pressure around the fluidic generator. The 
.cavity pressure p^ and ambient pressure p« were )nonitored. Figure 3 also 
Indicates the various electric outputs of the fluidic generator, such as the 
regulator output, the generator output, and the frequency output, which were 
recorded throughout this investigation. , 

Since the bench test setup, as shown in Figure 2, did not simulate the 
air flow around the generator, experiments were performed by placing the 
generator in the free- jet experimental test setup shown In Figure 4. The 
total stagnation inlet pressure, cavity pressure pc, and the pressure at the 
exit holes, along with the output of the fluidic generator were measured at 
various flight-simulated flow conditions. 

To determine the mechanism of operation of the fluidic generator, three 
1/8-in. fast-response pressure transducers with a rise time of 2psec and 


a frequency response frota 2 to 40»000 Hz were o>iip]oyed. One transducer was 
located at the base and the second was located flushed with the mouth of the 
tube to measure time varying pressure In the tube. The outputs of these 
transducers were cross-correlated on an all -digital -correlator. The third 
transducer was utilized to measure the cavity pressure fluctuations. Pigures 
5 and 6 show the physical locations of these transducers In the fluidic 
generator. 

To determine the dynamic response and the natural frequency of operation 
of the generator diaphragm-reed system* an accelerometer was mounted to the 
reed as shown In Figure 3. The output of this accelerometer was recorded on 
an rms meter and displayed on an oscilloscope. 



in. EXPERIMENTAL RESULTS 


3.1 Calculations of Laboratory Conditions to Simulate Free-Flight Conditions 
Given the altitude Z, the static conditions In front of the fluidic 
generator (denoted by subscript 1 in Figure 7) can be calculated using the 
U.S. Standard Atmosphere (Reference fl). To simplify calculations we use the 
geopotential altitude H, which equates the work required to lift a unit mass 
(under varying gravitational field) to Z with the work required to lift the 
same unit mass to H through a region with uniform gravity go “ 9.B0655 m/s2. 


Transformation: 

roZ 

From Z H H » — — ; ro “ 6,3S6,766m 

ro+Z 

- 0.0065H (KH <11, 000m 


T« 


216.65 

216.65 + 0.001 (H-20,000) 


-5.2559 

2B8.15 

/ 101,325.0 ( — ^ — ) 


11.000<H<20,000 

20,000<H<30,000 

0<H<1 1,000m 


( 1 ) 




1 22,632.1 
5,473.1 


22,632.1 exp [-1.5772 (11-11,000)10"'^] 

ll,000m<H<20,000m 

216.65 34.163 

(_- ) 2n.000m<H<30,000m 


( 2 ) 


Mo P« P«. 

3.4837 ~ 

T« 

where is the molecular weight and R'^ is the universal gas constant-^ In 

2 3 

these formulas T^ is obtained in K, p^ in newton per m , and p^ in kg/m . 

R* 

The sonic velocity at point (1) a« = x/yRT® » i"or air y “ ~ - 287.05 

Mo 


joule/ (kg K). 


- 5 - 


a„ «ylm.B7 T« « 20.046 Vt^, T„ In kelvins 


(3) 


Ml 


, U Is the fluidic generator velocity In in/s, 


Moving across the shock to point (2) 

P 2 2 y Ml - (y-1) ^ «1 -1 


Pi 


Y+1 


9 v/V"l Y-1 o 7 

P02 /(Y^l)Miy / Y+1 \ /<>Mi\ / 6 \ 

p 7 “(— ) i- ) '(-) (— ) 

2yMi-(y-1) 7Mi'-1 


1 

.y" 


7 

.2 ? 


5 

7 


2 (y-I)^! +2 Ml + 5 

M2 = ** 

2yMi - (y-1) 7Mi -1 


(4) 


T2 C2yMi - (y-1)3C(y-1) +2] (7M^-1)(M? + 5) 

^1 2 2 2 
(y+ 1) Ml 36Mi 

The important variable is the driving pressure Ap, which is the 
difference of the stagnation pressure Po 2 at the front of the fluidic 
generator minus the cavity pressure around the resonator (Ap ® Po2 " Pc)* To 
transform free flight to the bench test* it was assumed that Ap = Po2 - Pc ® 
P 3 - pe» (see Figure 3). Because of difficulties In estimating pc in free 
flight, a similar idea was tried by equating (Po - Ps) of po2 - Ps = P3 " P® 
(see Figures 3 and 7). 

Therefore, by knowing the altitude Z, and the fluidic generator velocity 
U, Ml = U/aco(Z) at that altitude can be calculated from relations (1) to (3) 
above. Using shock relation (4), the stagnation pressure in front of the 
generator po2 and therefore the laboratory conditions needed to simulate the 
flight can be computed. A computer program was developed, which Is attached 
in Appendix A along with some sample calculations. 
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3.2 Comparison of Bench Results With Free-Jet Results 

To determine the Influence of flight flow around the fluidic generator, 
bench test results (Figure 2) were compared with those obtained in the free- jet 
facility (Figure 4), The driving pressure (pq - p®) was kept constant. The 
pressure P 3 In Figure 3 Is assumed equal to the stagnation pressure Po in front 
of the generator in the free- jet facility shown in Figure 2. Typical results 
of such a comparison are shown In Figure 8. As Is evident In these results, 
for a given differential pressure Po - p„, the bench test output was always 
greater than that obtained by the free- jet test results. It should be noted 
that for this comparison the same generator was utilized. A comparison with 
two other generators was also made, and similar results were obtained. These 
results Indicated that the fluidic generator output did not depend upon the 
differential pressure po - p®. 

Upon careful analysis of other experimental data, I.e., cavity pressure Pc 
as a function of the driving pressure Po - p®, it was observed that for a given 
(Po “ P®)» cavity pressure p^ was always high in the free-jet case as 
compared with the bench tests. These results are shown in Figure 9. 

Therefore, in Figure 10, the fluidic generator output was plotted as a function 
0^ (Po " Pc) 3nd the agreement between free-jet and bench test results was 
excellent. 

From this study, it is therefore concluded that the most important 
parameter that controls the fluidic generator output is not Pq -P® but rather 
Po ~Pc* differential pressure between the total stagnation pressure in 
front of the generator and the cavity pressure. 


- 7 - 



3.3 Comparison of Laboratory Results With Flight Data 

Apart from getting the physical insight needed to model the fluid 
electric power generator from the above experiments, an attempt was made to 
correlate the laboratory tests with those obtained from flight data. The 
flight data was provided to JPL by Harry Diamond Laboratories. From the 
knowledge of the flight velocity and altitude, the laboratory flow conditions 
were calculated as discussed above In Section 3.1 and Appendix A. 

The results of such a comparison are shown in Figure 11. Since Pc was 
not available in the flight data, the static pressure behind the oblique shock 
was equated to pc. As is evident from Figure 11, the comparison with the 
bench test results is very poor. The discrepancy in such a comparison may 
have been caused by the assumption that cavity pressure pc in flight is equal 
to p„ around the generator close to the holes. Flight measurements of pc are 
needed to enhance our understanding to better correlate the flight data with 
laboratory test results. Furthermore, it should be noted that the absolute 
pressures in the bench facility and in flight are not identical. 
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3,4 Mass Flow Measurements 


To better understand the performance of the operation of the fluidic 
generator, detailed measurements of the mass flow through the fluidic 
generator were determined as a function of various operating pressures. The 
influence ov the exit hole area on generator performance was also determined. 

The calculations utilized to determine the mass flow rate through the 
fluidic device are shown in Appendix B. Figure 12 Indicates the mass flow 
rate normalized with cavity pressure and plotted as a function of pressure 
ratio Ps/Pc where pg is the static pressure at the holes, which under condi- 
tions of no external flow 1s the same as atmospheric pressure poa. Tables 1 - 4 
give the measured data along with the calculated mass flow rate and computed 
nozzle flow coefficients. The mass flow rate, and consequently, generator 
output drop for a given pressure ratio pg/Pc when the holes are closed 
partially. However, if the mass flow rate is normalized with stagnation 
pressure Po, then the results shown in Figure 12 collapse into a single curve 
as shown in Figure 13, Results of Figures 12 and 13 clearly indicate the 
importance of the cavity pressure p^ in controlling the mass flow rate through 
the generator. Any changes in the geometry of the holes will influence Pc 
which in turn modifies the mass flow rate and the output of the fluidic 
generator. 

To determine the influence of the ambient pressure pg = p® in the bench 
test on mass flow rate through the generator, the environment chamber 
shown in Figure 3 was employed. Figure 14 shows the typical results of such 
an investigation. For all the exit holes opened, the mass flow rate m was 
independent of the ambient pressure for a given stagnation pressure Po and Pc. 
These results strongly suggest that the mass flow rate under subsonic flow 
conditions through the generator critically depends upon the differential 
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pressure (Po - Pc). 

The influence of the Reynolds number calculated on the nozzle velocity 
coefficient is shown in Figure 15. These results showed that the changes in 
the exit hole area did not significantly modify the ram- jet flow responsible 
for generating pressure fluctuations in the table. Tables 1-4 give the 
output of the nozzle velocity coefficient for various operating conditions* 

3.5 Mechanism of Fluidic Generator Operation 

To correctly model the operation of the fluidic generator, important 
information was needed to determine the mechanism responsible for pressure 
fluctuations in the tube. As shown in Figures 5 and 6, the pressure 
transducer output was cross-correlated to determine any standing waves in the 
tube. 

The instant pressure transducer output at the mouth and base of the 
resonance tube of the fluidic generator is shown in Figure 16. To mount the 
transducer at the base, the diaphragm was replaced by a solid plate. Figure 
16 also shows the auto and cross-correlations of these microphone signals. As 
is evident, the peak of the cross-correlation signal was delayed by 1/4 the 
time period of the pressure fluctuation frequency in the tube. This frequency 
is the fundamental tube resonance frequency given by f=ia/4L, where L is the 
length of the tube. From these measurements, it is concluded that the 
pressure fluctuations, and consequently the oscillations of the diaphragm-reed 
system of the fluidic generator, result from an acoustic resonance phenomenon 
in the tube. These measurements did not support the Helmholtz type of 
oscillation phenomenon in the generator as has been commonly assumed in the 
past . 

The influence of mass flow rate on the output of the pressure 
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oscillations at the base of the resonance tube (so called since the frequency 
of pressure oscillations Is at the fundamental resonance tube frequency) was 
studied. Results In Figure 17 Indicated a linear relationship between the 
mass flow rate and the rms pressure fluctuations at the base of the cavity. 
This linear relationship suggests that the annular ram-air jet Is periodically 
swallowed by the tube at the frequency of resonance tube oscillation. This 
phenomenon, discussed in detail In Reference 9, Is called the regurgitant mode 
of resonance tube operation. Furthermore, since mass flow rate is 
proportional to (Po-Pc) ^s discussed in Section 3.4, the magnitude of the 
pressure fluctuations at the base of the resonance tube Is directly related to 
this differential pressure. 

3.6 Accelerometer Output 

Other experimental information needed to properly model the fluidic 
generator required us to determine the transient response of the fluidic 
diaphragm- reed subsystem of the fluidic generator. The accelerometer output 
was compared with the generator output at different pressures, as shown in 
Figure 18. Both outputs were at the resonance tube frequency. The transient 
response of the diaphragm-reed system is shown in Figure 19. The present 
measurements show the natural frequency of oscillation of the diaphragm-reed 
subsystem to be about 1412 Hz. The information shown in Figures 18 and 19 was 
utilized in an analytical model employed to compute the performance of this 
fluidic generator. 
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IV. ANALYTICAL RESULTS 


4.1 Governing Equations of Fluidic Generator Operation 

Since the operation of the fluidic generator performance critically 
depends upon the mass flow rate through the system, the governing equations 
which couple the annular ram-air jet, ring-tone jet-layer, cavity and 
resonance tube have been written in terms of nomenclature and Figure 20. The 
governing equations can be written as 



Annular Jet 


Jet Layer 
Resonance Tube 


Cavity 



, dPc 

where ( )' denotes acoustical oscillations and { — ) denotes the variation of 

dmj ^ 

the cavity pressure with the jet mass flow at the working (w) point of the 


fluidic generator. For detailed derivation of the equations, see Appendix C. 


For solutions of the form 



z » a+1u), in vector notation we write 


exp (zt) 


where 


a + 1w, M * 


2zL 

1 + R exp C- — ] 
a 


2zL 

1 - R exp C 3 

a 


One finds In matrix notation 



Introducing nondlmensional variables, 


At 

“J * Ij » PJ * 


Ay /dPc 


“H 


At 
Th- 
at 


/^\ 

W/ 


Irtij' w 


Ay/dPcN 

PH = -Vz’/ 
a dniH 


At 

“ TT — . PT * PT — 


w 


where 


2 L Vc 

C - — »K = — .PJ" 
a AyL 


mj 


my 


PH 


L - Length of the resonance tube 
Ay = Area of the resonance tube 
a = Speed of sound » average 
The governing equation can br written 
0 Kc -1 1 -1 

0 1 ajC+Pj 0 0 

0-1 0 apC+PH 0 

-1 1 0 0 ayC+Py, 

\ 1 0 0 0 M 


\ 


For solution, the determinant should be 
0 Kc -1 1 -1 

0 1 ayC+Pj 0 0 

0 -1 0 a^C+PH P “ 

-1 1 I 1 ayC+py 

1 0 0 0 M 


mn 

A 

my 


i.\ 

f 'Jc ' 
PJ 



zero 
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which results 1n the equation 

F(C)5 M + arC+PT K?(ajC+pj)(«HC+PH)+(«jC+PjH(aHC+PH)l 
+ («jC+Pj)(aHC+PH) ■ 0 , (B) 

l-*+Tlc(l+ — ) 


wHh H 


eC+Re"5 

eC.Re-C 


1+^+nC 


2 

paAp 


i'4) 




AjfD 


» h 




fi * - / — 

mo 


mo a 

fo L a 
where I) denotes the diaphragm-reed conditions. 

For the Solution of the above equation, we need to calculate the 


constants first. 

4.2 Calculations of Constants 

For these calculations, the dimensions of the fluidic generator have been 
used as follows; 

Diaphragm-Reed 

For these calculations, data in Figure 19 has been utilized, 
mo a 

T, = « 22.5 

L 


r = 0.946 
«}i 1s unknown 
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Jet 


Holes 


or 


h 

a, * * .5243 

^ Aj L 


/Pc \ 

7CH- ) ] 

n pc \po / 

■ .006 — — u—- *_ 

m /pc\2/7 
5-6 (-) 

'Po' 

" - for Pc/Po < 0.528 

At 


for p^,/p^ 


«K 


.0936 


AhL 


2 / 7 . 


pc 7 Cl-(Ps/Pc) 3 Pc 

PH * *006 — for ~ 

Po 




2 -(Ps/Pc) 


a Pc Ps 

Ph = 0.0006 — for — < 0.528 

m Pc 


Jet Layer 


hAy 

OCT * = 0.0132 

AWL 


PT = 


Ay ifi 


p a A^ w 2 


349 m 


Cavity 


K - = 14.5831 

AtL 

Independent Variables 


C = (a+iw) I = 1.095 X 10 ’^ (cr+lw) 
a 


> 0.528 


> 0.528 



4.3 Numerical Results 


A computer program was developed to solve the equation (5) as developed 
In 3,1 above with the above physical constant. The computer program Is listed 
In Appendix D. The Input variables were the working conditions as follows: 

Po stagnation pressure In front of the oscillations 
Pc cavity pressure 

Ps static pressure close to the exit holes 
m mass flow rate 

The solution of equation (5), F(c) = 0 was found for the following cases; 

(a) Fixed resonant frequency of the diaphragm (q) for various 4> at 
the operating conditions of the fluidic generator. 

(b) Fixed (j) and varied r for these three operating conditions. 


The output of these results are shown In the following table. 






1 



o)L 


0 L 

Po 

Pc 

Ps 

m 

n 

Q 


a 


a 







41=^0 

4.=1 

p=o 

p=l 





I . I j 



1.63 

wmum 


20.38 

15.147 

14.32 

.0270 

22.527 

0.946 

4.67 

4.68 

-.093 

-.093 







1.59 

1.70 

-.1^7 






22.527 

1.5 

4.61 

4.68 

-.138 

-.137 







1.69 

1755“ 

-.131 

-.133 





22.527 

.946 

4.67 

4.68 

-.138 

-.137 







1.59 

isai 



26.559 

16.093 

14.32 

.0397 

22.527 

.5 

4.61 

4.67 

-.138 

1 

a-* 

CO 





■■1 

mn 

in 

1.68 

“Tnr 

-.137 










“ -274’4’ ' 





1 

.946 

4.61 

4.87 

-.130 

-.172 







1.60 

T.63 

-.175 

-• 169 

32.781 

17.082 

14.32 

.0505 

22.527 

.946 

4.66 

4.67 

-.175 

-.175 
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It is evident from the numerical results that 

(a) The natural frequency of the system is controlled hy resonance 
conditions of the tube I.e,, wl/a « n/2, 3^/2 

(b) Oscillations are always damped 

(c) To keep the system in operation, forcing of the systeni is 
needed. 
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V. DISCUSSION AND CONCLUSIONS 


The experimental work performed to develop the mathematical model to 
predict generator performance has shov^n that the mode of cavity operation is 
governed by the standing acoustic resonance phenomenon in the tube. In this 
mode of operation, the annular jet is periodically swallowed and subsequently 
discharged by the resonance tube at the fundamental tube resonance frequency 
(Reference 9). The basic mechanism of conversion of flow energy to acoustic 
power (and subsequently to mechanical motion) by the fluidic generator is 
critically dependent on the jet mass flow and the geometry of the gap between 
the nozzle exit and the resonating cavity. The basic understanding of the 
physical processes in this region under a simulated high altitude flow 
configuration is essential for designing the fluidic generator with optimum 
performance. More work is needed in this area to improve the efficiency of 
the conversion of fluid energy to electric power under high-altitude flight 
conditions. 

The comparison of laboratory test results with the fligr.t test data 
(Figure U) shows that the existing laboratory tests do not correctly simulate 
the performance of the fluidic generator at altitude, A close look at the 
flow field in the annular jet showed that the flow becomes supersonic inside 
the generator in flight,, as compared with subsonic flow in the existing 
laboratory simulation facilities. To correctly simulate flight, a vacuum 
facility is needed to simulate both the high-altitude mass flow through the 
fluidic electric generator and the static flow pressures in and around the 
fluidic generator. 
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APPENDIX A 

CALCULATIONS OF LABORATORY CONDITIONS TO 
SIMULATE FLIGHT CONDITIONS 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


fUkmAN 

AIHl.FTN 


/ 

/ 

,/ 


0001 

3002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 


0013 

0014 

0015 


0016 

0017 

0018 

0019 

0020 
0021 
0022 
0023 


0024 

0025 


002-51 12 : 24:12 15~0CT*-79 PAGE 1 

7TRi!BL0CKS/WR 

C THIS PROGR^AM DETERMINES CONDITIONS AT THE END OP THE 
C FLUIDIC GEN. GIVEN THE VELOCITY AND ALTITUDE OF THE MISSILE, 

C FROM THESE CONDITIONS AT THE END OF THE FLUIDIC GEN, 

C THE CONDITIONS UPSTREAM IN THE AIR STRAIGHTENING SECTION ARE DETERMINED 
C PS = STATIC SURFACE PRESSURE OF OGIVE AT PT, 2, 

CP PRESSURE <N/M»*2> AT DIFFERENT POINTS. 

C PO = STAGNATION PRESSURE <N/M*»2) AT DIFFERENT POINTS. 

C T = TEMPERATURE < KELVIN) AT DIFFERENT POINTS, 

C TO a TOTAL TEMPERATURE <KELVIN) AT DIFFERENT POINTS, 

CM a MACH NO. AT DIFFERENT POINTS. 

C PSI a P IN (PSD 

C TR a T IN <R> 

C PT. 1 a CONDITIONS UPSTREAM OF BOU SHOCK WAVE. 

C, 2a CONDITIONS DOWNSTREAM OF BOU SHOCK (ACTUAL FLIGHT COND.) 

C JUST UPSTREAM OF FLUIDIC GEN. 

C PT. 3 a CONDITIONS JUST UPSTREAM OF FLUIDIC GEN. IN LAB. 

C APPARATUS. 

C PT. 4 a CONDITIONS IN THE AIR STRAIGHTENING SECTION OF LAB, 

C APPARATUS. 

DIMENSION HEAD(2) 

DOUBLE PRECISION P(4> »T(4) »PSI(4) »TR(4> rP0(4) »T0(4) » 

1 CP(2)rCV(2>fM(4)rK»A3»A4»RH0(4)»PS»X 

DATA A3/1.266769/»A4/182. 41469/ 

DATA HEAD/4HAIR r4HHEL./»CP/3.5f 2.5/»CV/2.5»l .5/ 

DATA XSAVE/1 . OE+20/ r PCONV/1 , 45038E-04/ 

TYPE 'INPUT VELOCIYT(M/S) rALTTTUDE(M)' 

ACCEPT *»V»Z 

CALL ATM76(Z»T(l)rP(l>»RH0(l>) 

A a 20.046803 * SQRT(T(1)) 

M(l) a V / A 

PO(l) a P(l) * <1. + .2*M(1)»M(1)>**3.5 
rod) a T(l) * <1. + .2*M(1)*M<1>) 

C 

C CALC. MACH NO. OF MISSILE. 

C 

A a 20.046803 » SaRT(T(i>> 

V a M(l) ♦ A 
IF (M(1>.6T.1.)G0 TO S 
C 

C SUBSONIC FLIGHT 

C 

T(2) a T(l> 

V(2) a p(i) 

M(2) a M(l) 

P0(2) a P0(1> 

T0(2> a T0(1) 

PS a P(2) + .23 * (P0(2) - P(2)> 

GO TO 6 

5 CONTINUE 

C 

C CALC. OF CONDITIONS AT FLUIDIC GEN, FROM ATMOSPHERIC CONDITIONS 
C AND VELOCITY OF MISSILErSAME PROBLEM AS CALC. DOWNSTREAM CONDITIONS 
n t 'F A NORMAJ FHOCK WAVE FROM .UPSTREAM COND I T IIINS . 

X a M(l) * M(i) 

P(2) a p(i) « <7.*X “ 1. ) / 6. 
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FORTRAN IV- 
ATNl^FTN 

0026 

0027 

002B 

0029 


PLUS V02-51 

/TR: BLOCKS/UR 


12i24il2 


15-0CT-79 


\ 

PAGE 2 


0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 


0046 

0047 


P0<2^s P(l) $ <1.2*X)**3.5 * <P(1)/P(2))*#2.5 

MC2) = SQRT<(X+5.) / <7.*X-1.)> \. 

T<2) » T<1> * <7,»X-1») * <X+5.) / <36, *X) 

T0(2) = T<2) $ <1, + .2*M<2)*H<2)) 

CALC OF CONDITIONS ^UST UPSTREAH OF FLUIDIC GEN, FOR 
THE LAB APPARATUS FROM THE ACTUAL UPSTREAM CONDITIONS 
REQUIRED CONSTRAINTS ARE {DRIVING PRESSURE OF THE GEN. 

MUST BE THE SAME. 

{THE STAGNATION TEMP. OF THE ATMOSPHERE 
JUST BEFORE THE FLUIDIC GEN, <PT. 2) MUST EQUAL THE STAGNATION 
TEMP. At PT, 3. 

DO 40 Leif 2 
K = CP(L> / CV<L> 

P<3) = 101325.0 

P0<3) = P0<2) + P<3) - PS 

M<3> = SQRT< 2./(K - 1.) $ < (P0<3)/P<3> )**< <K - 1.>/K> - 1.) ! 
STORE = ,5*<K - l.>#M<3MtM<3) 

TOO) ® T0<2) 

TO) “ T0<3> / <1. + STORE) 

EXPON « .5 * <K + 1, > / <K - 1.) 

CRIT » A47(A3»M(3)) « (1. STORE)/( .5«<K-M . ) )*«EXPON 
X » (1. / (.5»(K+l.)>)i»EXP0N 
M<4) = X / CRIT 
P0(4> » P0(3) 

T0<4> = T0(3) 

P<4) = P0<4) / <1, + ,5*<K-1.)»M<4)*M<4))**<K/<K-1)> 

T<4) != T0<4) / <1. + .5«(K-1.)»M<4)«M<4>) 

PRINT OUT RESULTS. 

TYPE 100fHEAD<L)»VfZ 

) F0RMAT</»1X» 'GAS IN LAB APPARATUS { ' »A4/2Xr 'VELOCITY = 'iF8.2i 


' <M/S) '/2X» 'ALTITUDE = '»F9. 
31X»'STAG, STAG 
IXf'PT. TEMP.<K) P <N/M**2) 
' P <PSI) MACH NO,') 


»F9.2»' (M)'f// 
STAG.'f/ 

**2) TEMP.<R) 


TEMP,<R) 


(PSD 


0048 


DO 30 1=1 r4 

0049 


TR<I) = 1,8 * T(I) 

0050 


PSKI) = PCONV ♦ P(I) 

0051 


TRO = 1.8 * T0<I) 

0052 


PSIO = PCONV ♦ PO(I) 

0053 

30 

MRITE<6»110)IfT<I) »P<I)rTROrTR(I)fPSIOfPSI(I)fM<I) 

0054 

110 

F0RMAT<lXrI2f2X»F8.2»F13,2f2F10.2.2F11.5»F12,6) 

0055 


X = <P0(2) - PS) » PCONV 

0056 


URITE(6f 120)X 

0057 

120 

F0RMAT<2X» 'DRIVING PRESS. <P(4)-P<ATM. ) > = 'fF7,3) 

0058 

40 

CONTINUE 

0059 


STOP 

0060 


END 
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12J28J52 


15-0CT-7V 


PAGE A 


FORTPAN IV-PI.US V02-51 
A TM 1 ♦ FTN / TR { BLOCKS/ WR 

0001 SUBROUTINE ATM7<j»<ZrT»P»RH0) 

C 

C SUBROUTINE ATM76 GENERATES THE ATHOSPHERIC CONDITIONS 
C TrP»RHO GIVEN THE ALTITUDE ACCORDING TO THE U.S. STANDARD 
C ATM0SPHERErl97A. 

C 


0002 


DOUBLE PRECISION TrPt RHO,MWTf G»RrROrZ»H 

0003 


DATA R0/63S6766 ♦ / » NWT/28 . 9655/ 1 Q/9 . 80655/ » R/8314 . 32/ 

OOOA 


H « RO ♦ Z / <R0 + Z) 

0005 


IF <H.LE,11000«)G0 TO 10 

0006 


IF < H.LE. 20000 * )60 TO 20 

0007 


IF ( H.LE. 32000. ) GO TO 30 

0008 


TYPE *» 'ALTITUDE IS TOO HIGH' 

0009 


RETURN 

0010 

10 

T s 288,15 - .0065 » H 

0011 


P “ 101325.0 * (288.15 / T)*»< <~1 . )*G*MWT/<R*.006S> ) 

0012 


GO TO 40 

0013 

20 

T = 216,65 

0014 


P « 22632.4 t EXP<(-1.)»G*MWT*<H-11000,)/<R*216,65)> 

0015 


GO TO 40 

0016 

30 

T s 216.65 -f .001 « (H - 20000.) 

0017 


P = 5475.05 ♦ <216.65 / T)**<G*HWT/R) 

0018 

40 

RHO = HUT*P / <R*T) 

0019 


RETURN 

0020 


END 


FORTRAN IV-PLUS V02-51 12:28J52 15-0CT-79 PAGE 5 

ATHl , FTN /TR J BLOCKS/WR 

PROGRAH SECTIONS 


NUHBER 

NAME 

SIZE 



ATTRIBUTES 




1 

$C0DE1 000474 

158 


RUfIrCONfLCL 




2 

iPDATA 000136 

47 


RUrDfCONfLCL 




4 

4VARS 

000050 

20 


RUrDfCONrLCL 




ENTRY 1 

POINTS 








NAME 

TYPE 

ADDRESS 

NAME 

TYPE 

ADDRESS NAME 

TYPE 

ADDRESS 

name 

ATM76 


1-000000 







VARIABLES 








NAME 

type 

ADDRESS 

NAME 

TYPE 

ADDRESS NAME 

TYPE 

ADDRESS 

NAME 

C 

R«8 

4-000010 

H 

R*8 

4-000040 MUT 

R»8 

4-000000 

P 

RHO 

R«B 

P-000010* 

RO 

R«8 

4-000030 T 

R*a 

F-000004* 

Z 
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KUN ATMl 

INPUr V£LOCIYT(M/S) »ALTITUDE<H) 
I $'9. 14209. 


GAS IN LAB apparatus: AIP 
VELOCITY » 199.00 <H/S) 

ALTITUDE » 4209.00 (H) 





STAG. 


STAG. 



PT. 

TEHP.CK) 

P <N/H»*2) 

TENP.<R> 

TEHP. <R) 

P <PSI) 

P (PSD 

HACH NO. 

1 

240.81 

40000.78 

504.93 

469.46 

11.22968 

8.70239 

0.614676 

2 

240.81 

60000.78 

504.93 

469.46 

11.22968 

8.70239 

0.614676 

3 

270.73 

101325.00 

504.93 

487.31 

16.64199 

14.69597 

0.425257 

4 

280.52 

114741.60 

504.93 

504.93 

16.64199 

16.64189 

0.002850 

DRIVING PRESS 

<p{4)~P(ATH 

.)) => 1 

.946 




GAS 

IN LAB APPAimUS: HEL, 






VELOCITY » 

199.00 <H/S> 






ALTITUDE » 

4209,00 <H) 









STAG, 


STAG. 



PT. 

TEHP,(K> 

P <N/H**2> 

TEMP.<R) 

TEHP. (R> 

P (PSD 

P (PSD 

HACH NO. 

1 

260,81 

60000,78 

504.93 

469.46 

11.22968 

8.70239 

0.614676 

2 

260.81 

60000.78 

504.93 

469.46 

11.22968 

8.70239 

0.614676 

3 

266.91 

101325.00 

504.93 

480.43 

16.64199 

14.69597 

0.391152 

4 

280.52 

114741.62 

504.93 

504.93 

16.64199 

16.6418? 

0.002585 


DRIVING PRESS. <P(4)-P< ATH. ) ) » 1.946 

TTO — STOP 

>RUN ATHl 


INPUT VELOCIYT(H/S>rALTITUDE<H) 
176.9.4150.5 


GAS IN LAB apparatus: AIR 


VELOCITY = 

176*90 (H/S> 






ALTITUDE = 

4150.50 (H> 









STAG. 


STAG. 



PT. 

TEHP.(K) 

P (N/H**2) 

TEHP.(R> 

TEHP.(R) 

P (PSD 

P (PSD 

HACH NO. 

1 

261.19 

60461.39 

498,17 

470.14 

10.73977 

8.76920 

0.546016 

2 

261.19 

60461.39 

498.17 

470.14 

10.73977 

8.76920 

0.546016 

3 

269.10 

101325.00 

498.17 

484.38 

16.21331 

14.69597 

0,377305 

4 

276.76 

111786.15 

498.17 

498.17 

16.21331 

16.21324 

0.002548 

DRIVING PRESS 

, (P(4)-P(ATH. 

>> - 1 

.517 





GAS IN LAB APPARATUS: HEL. 
VELOCITY - 176.90 (M/S> 

ALTITUDE « 4150.50 <M) 


PT. 

TEMP. (K> 

P (N/H**2) 

STAG. 

TEHP.(R) 

TEHP. (R) 

STAG. 

P (PSD 

P (PSD 

1 

261,19 

60461.39 

498.17 

470.14 

10,73977 

8,76920 

2 

261.19 

60461.39 

498.17 

470.14 

10.73977 

8.76920 

3 

266.10 

101325.00 

498.17 

478.97 

16.21331 

14.69597 

4 

276.76 

111786.16 

498.17 

498.17 

16.21331 

16.21324 


driving press. <P<4)-P<ATH. >) = 1.517 

TTO — STOP 


HACH NO. 
0.546016 
0.546016 
0.346783 
0.002315 
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RUN ATHl 

INPUT VELOCiyT<M/S) » ALTITUDE <H> 
18?if7»319S*4 

GAS IN LAB apparatus: AIR 
VELOCITY « 185.70 <M/S> 

ALTITODie = 3195.40 <N) 




STAG. 


STAG, 



PT. TEHP.IK) 

P 

TEHP.<R> 

TEMP. <R) 

P <PSI> 

P (PSD 

MACH NO. 

1 267.39 

68397.19 

512,19 

461.30 

12.33323 

9.92019 

0.566492 

2 267.39 

68397.19 

512.19 

481.30 

12.33323 

9,92019 

0.566492 

3 275.04 

101325.00 

512.19 

495.06 

16.55401 

14.69597 

0.415938 

4 284.55 

114135.07 

512.19 

512.19 

16.55401 

16.55392 

0.002792 

DRIVING PRESS 

1. <P<4)-P<ATH. 

)) * 1 

.858 




GAS IN LAB apparatus: HEL. 






VELOCITT * 

185.70 (N/S> 






ALTITUDE = 

3195^40 (M) 

STAG. 


STAG. 



PT. TEHP.(K) 

P (N/M**2) 

TEMP, (R) 

TEMP.(R) 

P <PSI) 

P (PSD 

MACH NO. 

1 267.39 

68397.19 

512.19 

.481.30 

12.33323 

9.92019 

0.566492 

2 267.39 

68397.19 

512.19 

481.30 

12.33323 

9,92019 

0.566492 

3 271.32 

101325.00 

512.19 

488.37 

16.55401 

14.69597 

0.382521 

4 284.55 

114135.08 

512.19 

512.19 

16.55401 

16.55392 

0.002533 

DRIVING PRESS 
TTO — STOP 

. <P(4)-P<ATH. 

>) '- 1 

.858 




>RUN ATMl 







INPUT VELOCIYT<M/S>rALTITUDE<H> 





203.4rl615.4 







GAS IN LAB apparatus: AIR 






VELOCITY -- 

203.40 (H/S) 






ALTITUDE - 

1615.40 (M) 

stag. 


STAG. 



PT. TEMP.<K) 

P <N/H**2) 

TEMP.<R> 

TEMP. <R) 

P <PSI) 

P (PSD 

MACH NO. 

1 277.65 

83369.21 

536.84 

499.77 

15.53180 

12.09170 

0.608913 

2 277.65 

83369.21 

536.84 

499.77 

15.53180 

12.09170 

0.608913 

3 284.45 

101325.00 

536.84 

512.01 

17.34485 

14.69597 

0.492382 

4 298 , 24 

119587,41 

536.84 

536.83 

17.34485 

17.34472 

0.003261 

DRIVING PRESS 

. <P<4)-P<ATM. 

)) = 2.649 




GAS IN LAB apparatus: HEL. 






VELOCITY - 

203.40 <M/S) 






ALTITUDE = 

1615.40 <H) 

STAG. 


STAG. 



PT. TEHP.<K) 

P <N/H»*2) 

TEMP. <R) 

fEMP. <R) 

P (PSD 

P (PSD 

MACH NO. 

1 277.65 

83369.21 

536.84 

499.77 

15.531S0 

12.09170 

0.608913 

2 277.65 

83369.21 

536.84 

499.77 

15.53180 

12.09170 

0.608913 

3 279.11 

101325.00 

536.84 

502.40 

17.34485 

14.69597 

0.453438 

4 298.24 

119587.43 

536.84 

536 . 83 

17.34485 

17.34472 

0.002947 

DRIVING PRESS 
TTO — STOP 

* (P<4)“P<ATM. 

>> = 2 

.649 






APPENDIX B 


CALCULATIONS OF MASS FLOW 
THROUGH THE FLUIDIC GENERATOR 


i 

i 

i 

i 

1 

i 


i 


I. H^^SS Flow CALCULATION THROUGH THE FLUIDIC GENERATOR 

It 1$ Uesirfible to know the relationship between the absolute total 
pressure Po» the absolute cavity pressure Pc» and the absolute skin pressure 
Just outside the exit hole p^. To obtain these relationships, assume Ideal 
Isentropic flow of air through a nozzle and multiply the mass flow rate 
expression by Cy, the velocity coefficient: 

3 


. , , Jzy PlPl 

m » CyA2\“* - — ■ . " ■ 


If State (0) Is the stagnation condition »then for A2«Ai, the mass flow 


rate simplifies for air to 



Assuming the Ideal gas relationship ° Pr> RTn where R Is the gas constant, 

. _ ho! /P >.2/^ 


RTq 'Po^ \On/ 


m « (2.0546 CyApo//To 

where A " In^, Pq ~ Psta, T " °R, m " Ibm/sec. 

Cy Is usually assumed a function of the Reynolds number of the orifice 

or nozzle. This Implies that the mass flow rate from a high pressure to a low 

pressure through some arbitrary small opening Is a function of the dlmenslon- 

md 

less quantities y » ratio of specific heats; R^ = — ; p/p©. The mass flow 
rate will also depend on the stagnation pressure and temperature, Pq and Tq. 
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The converging section of the bench setup allows for a smooth, gradual 
convergence of the air just before It flows Into the generator. The mass flow 
rate was calculated from the pressure ratio relative to this section, assuming 
that flow through this Section Is isentropic. This approximation should give 
an error In the range of where the positive error indicates that the mass 
flow rates calculated are greater than the actual mass flow rates. A simple 
procedure to get the actual flow rate Is to calibrate our setup using a 
flowmeter and obtaining the relationship of Cy vs Reynolds number. The 
equation presently used to calculate the mass flow rate Is 



li. VELOCITY COEFFICIENT CALCULATIONS 


There are two velocity coefficients: Cyj, the expansion from the ram 
pressure to the cavity, and Cy^^, the expansion from the cavity through the 
exit holes to the atmosphere* Each of these coefficients Is a measure of the 
effective area that Is reduced by the Influence of the boundary layer at the 
edge of each orifice. 

In the present Investigation, Cy was studied as a function of Re only. 
The equation used to calculate C^j for the expansion from po-^Pc Is 

m Vfo 

CyJ - 

Similarly, 


mVTo 
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The variables used to calculate and CyH are described and their 


origins are noted 

In the table below. 


k il fn il miHH 

Description 


PoiTo 

Stagnation conditions 

Experiment measurement 

# 

m 

Hass flow rate 

Calculated from smooth 
converging section 

Aj.Ah 

Cross-sectional areas 

Drawings 

Ps 

Skin pressure 

Experiment measurement 

Pc 

Cavity pressure 

Experiment measurement 



APPENDIX C 


DERIVATION OF THE EQUATIONS OF 
FLUIDIC GENERATOR OPERATION 


DERIVATION OF EQUATIONS OF FLUID GENERATOR OPERATION 


The operation of the fluidic generator has been divided into five 
components consisting of the annular jet, the jet layer, the resonance tube, 
the cavity, and the holes (see Figures 1 and 7). A linear differential 
equation has been written for each of these five components that relate the 
mass flow and pressure fluctuations In different parts of the generator. 

Because of the linearity of the equation we look for solution of the form 

A 

p' = p exp(a + 1(o)t 

A 

»fi‘ = m exp(a + 1o))t (C.l) 

A 

where w is the frequency of oscillations, cr the damping coefficient and p and 

A 

m are the amplitude of the oscillations, A positive value of a indicates 
exponential growth at a particular frequency w. From the homogeneous system 
of equations, the amplitude ratios can be obtained. These ratios contain a 
phase factor giving the phase differences between the various pressures and 
mass flow through the fluidic generator. 

The outcome of the system of these equations Is a function a = a(w) which 
gives the damping or growth of acoustical waves of a given frequency and 
amplitude, and the phase relationship between pressures and mass flow at 
various locations, 

A, Annular Jet 

The annular jet Is a convergent nozzle; the relation between mass flow 
and pressure ratio (total pressure before the contraction divided by the 
static pressure at the exit) is well known. In this analysis we have assumed 
that the exit pressure equals the cavity pressure, thus 
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where Po» Pc are the total and cavity pressures respectively. The second 
equation corresponds to shock conditions at the exit. The mass flow 
coefficient Cyj incorporates possible viscous effects. Typical values of 
these coefficients are given in Table 1 to 4. 

Acoustical fluctuations in the cavity induce fluctuations of the mass 
flow which are given by the relationship. 

dmj 
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1 P’ 
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where ( )' denotes acoustical fluctuations and the derivative ^ is to be 
evaluated at the operating point 
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The previous relation does not take into consideration the Inertial 
effects on the relationship between the acoustical pressure and mass flow. In 
order to Incorporate these effects, an Inertial term was Incorporated into 
the equations. It Is assumed that the fluctuating pressure 1s In part utilized 
to accelerate the mass flow through the nozzle. The resulting equation is 



From dimensional considerations Ij will be a function of the operating 
conditions and the geometry of the jet nozzle. 

h 

~ fj geometry) 

Aj 

where h Is the width of the ring jet, Aj its area and fj varies from a value 

of order one at Mj = 0 to zero at Mj = 1. For the purpose of this analysis 

1 for Mj < 1 
fj(Mj) - _ 

0 for Mj > 1 


Seeking a solution as Indicated in equation (C. 1), the equation for 
annular jet mass flow can be written as 

/dpc\ - 

[Ij (cj+iw) + ( — ) ] mj = Pc 
Vdifij/ 

w 

B. det Layer 

The jet layer Is the region of interaction between the entering ram- jet 
and the lip of the resonant tube. In first approximation, the pressure 
inside and outside the ring jet layer can be considered uniform since the 
acoustical wavelength is much longer than the distance between the jet exit 
and the tube lip. The pressure Inside (pj) and outside (pc) need not be the 
same Instantaneously. Phase differences In their oscillations, however, will 

- 35 - 


result in a deflection force on the jet layer which will cause a radial 
displacement and therefore an oscillating mass flow into the tube (mj). 

Detailed analysis of this interaction is an extremely complex problem and 
is beyond the scope of current research. A phenomenological model was therefore 
developed in which the inertia and the curvature effects of the jet layer have 
been incorporated. The present model assumes curved streamlines with constant 
radius of curvature as sketched. Furthermore, the jet width is assumed small 
compared with the radius of curvature of the layer and the ring jet. 



X U 


Momentum conservation requires that the pressure differential be balanced 
by the inertia and the centrifugal force. The inertia term is given by 



where 


X 

V = q - 
Rj 


we find 

Ij 


where h is the jet width, W 


I 



q 


0 


X 


a 1 d(l/Rj) 

— f— )dx = phWq — 

at Rj' 2 dt 


its perimeter, and q the jet exit velocity assumed 
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constant and Independent of x. The centrifugal term is given by 



and finally the pressure force Is (Pj - Pc) AM, 


Thus A ^ 

d(- ) 

A 'Rj' q2 

P h q _ + p h =; (pj _ p ) 

2 dt Rj 

The radius of curvature is related with the tube mass flow m‘ by the equation 

ffi' J 

2 T 

1/Rj = . — 

A^ pqW 


The relation between the tube mass flow, the pressure in the cavity and the 
tube Is 

h d rfit 2(m)w 

— = P' - P' 

AW dt T c T 

where {m)yy is the mass flow through the system at the operating conditions. 
For sinusoidal solution we find: 



C. Resonant Tube 

The most important aspect of the operation of the fluidic generator is 
the energy transfer in the resonant tube. An adequate model of the system has 
to incorporate its two functions as a) the resonant tube system and b) the 
"acoustical transmission line" phenomenon through which the energy is 
transferred to the diaphragm. Depending upon the acoustical wavelength, a 
variety of tube resonant modes can be excited. At low frequencies, i.e, the 
wavelength is much longer than the diameter of the tube, only the longitudinal 
modes are excited. These modes are characterized by plane waves moving back 
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and forth through the tube. The difference in energy carried by the waves 
travelling in one direction compared to those moving in the opposite direction 
is the amount of energy transferred through the tube, The resonant 
characteristics of the tube are determined by the matching conditions at both 
ends which in turn determine the amount of energy transmitted. 

In the present model, the relationship between the pressure and the mass 
flow at the entrance of the tube is a function of the amount of energy 
transferred through the tube. 



We consider the superposition of right and left propagating acoustical waves. 
Their amplitude and phase at the diaphragm end are related by the reflexion 
coefficient R(w) which in general is a complex function of the frequency. 

In general, if the pressure field associated with the right propagating 


then at any x 


Pr - PR(t) 


PR(t.x) » Pr (t-— ) 
a 

where a is the speed of sound. Similarly for the left propagating wave 

X 

PL(tpX) = PL(t+-) 
a 

At the diaphragm end, the reflexion coefficient R is defined as 

R , »-/a) 

PR(t- L/a) 

For wave forms given by 

A 

PR(t) = Pr exp[(a+iw)t] 

Pl(t) = Pl exp[((j+iu))t3 

we find 


R = — exp [2(o+iu)) L/a] 

pR 


and consequently at x=o the pressure fluctuation is given by 

A A A 

Pt = (Pr + P|_) exp[(o+iw)t] = p exp[(o+iw)t] [1+R exp {-2(ar+io))L/ao}] 

R 

Associated with the right and left propagating waves there is an acoustical 
velocity field given by 


Ur = Pr/p3 
Ul = - Pl/p«> 

Therefore, the velocity at the entrance of the tube is given by 

.A 

Pr 

Uy = — - exp[(as'ito)t] {1-R exp(-2(a+ia)) L/a)} 

Pa 

The mass flow m' = p Uj Aj is 

^ A 

m' = my exp[(a+iti))t] 

T fs 

Pr 
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The relationship between my and pj Is 

•"T At 1 - R exp C-2(cy+1w) L/a] 


PT a 1 + R exp [-2(o+l4») L/a] 

The reflection coefficient R Is determined by the Interaction of the 
acoustical field with the diaphragm. The diaphragm can be modelled as a 
second order mechanical system. If xq Is its position then Its notion Is 

d?XD dxg 

described by the equation P*n An ® hIq — + f p — + kpxo 

dt2 dt 



where mp and kp are the mass and spring constant of the diaphragm-reed 

combination, and fp is the friction coefficient which also incorporates the 

effects of the energy transmitted into the electrical part of the system. The 

dXD 

fluid velocity at the diaphragm is -—which results In a velocity at the end 

dt 

of the tube U', 

t Ap Ap dxp 

U' , — U' = — — 

E At 0 At dt 

Since p'e at the end of the tube equals p^, we find 


= At/Ad 


The tube equation is given by 
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E UR 
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dxp 

— + kpxp 
dt 


- 40 - 


The reflection coefficient can be written as; 
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( 1 - 
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Arfo 


dXD 

) — + mo 

dt 
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(1+ ) + mo 

Ayfo dt 

For sinusoidal oscillations, one can write 
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where z = a+lw 

Introducing the nondimenslonal parameters 
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The term coefficient can be written as 
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+ kpxp 


+ kpxp 


D. Cavity 

Both the jet and the tube are located Inside a cavity which plays a role 
In the acoustical properties of the system. For the frequencies of fluidic 
generator operation, this cavity Is operating as an accumulator because the 
wavelength Is very large compared with the characteristic dimension of the 
cavity. The acoustical pressure Inside the cavity can also be considered 
uniform. Furthermore, the characteristic time of the oscillations is short 
enough to neglect the heat transfer from the fluid to the wall. Therefore the 
cavity can be modeled as an Isentropic accumulator for which the relation 
between the mass flow and cavity pressure Is given by 

dPc yRTc 

where y 1s the ratio of specific heats, R the gas constant, Tc the gas 
temperature and Vq the cavity volume. 

This equation gives for oscillatory waves 

yRTc , 

(a+1w) Pc = - — % 

Vc 


E. Holes 


The equation modeling the behavior of the exit holes Is Identical to that 
governing the behavior of the annular jet. It should be noted that for the 
exit holes, the acoustical fluctuations depend on the total pressure upstream 
of the orifice rather than at the exit pressure; thus; 
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Therefore the linearized equation Is 
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Ih = — "fH (Mex» geometry); fH (Mgx. geometry) = 1 

iidH 

fH has been assumed independent of the exit Mach number. 

Thus for sinusoidal oscillations, we find 
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IH (a+iw) + liiH “ Pc 

dmH • 
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APPENDIX 0 


COMPUTER PROGRAM FOR MODELING THE 
FLUIDIC POWER GENERATOR 
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Table 3. Data for 3/4 Exit Holes, Open Area 
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Table 4. Data for 1/2 Exit Holes, Open Area 
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OUTPUT WIRES 
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Figure 1. Fluidic Electric Power Generator 
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Figure 3. Instrumentation of Fluidic Electric Power Generator 
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Figure 4. Free-Jet Test Setup 
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Figure 6. Pressure Transducer Location at the Base and the Mouth of the 
Resonance Tube 
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BOW SHOCK 



Figure 7. Calculation of Bench Flow Conditions 
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Figure 8. Comparison of Laboratry Measurements with Flight Data 
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Figure 9. Cavity Pressure vs Free-Stream Dynamic Pressure 
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Figure 10. Comparison of Free-Jet Test Data with Bench Test Results 
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Figure 11. Comparison of Laboratory Measurements with Flight Data 
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Figure 13. 


Mass Flow as a Function of Pj./Pq Various Exit Areas 
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Figure lb. Velocity Coefficient vs Reynolds Number for Various Exit Areas 
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Figure 16. Pressure Transducer Output 
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Figure 17. RMS Output of Pressure Transducer as a Function of Mass Flow Rate 
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Figure 18. Oscilloscope Traces of Fluidic Generator 
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Figure 19. Transient Response of Diaphragm-Reed Subsystem 
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Figure 20. Modeling of Fluidic Electric Generator 








